Structure and stability of self-assembled actin— lysozyme complexes in salty water 
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Interactions between actin, an anionic polyelectrolyte, and lysozyme, a cationic globular protein, 
have been examined using a combination of synchrotron small-angle x-ray scattering and molec- 
ular dynamics simulations. Lysozyme initially bridges pairs of actin filaments, which relax into 
hexagonally-coordinated columnar complexes comprised of actin held together by incommensurate 
one-dimensional close-packed arrays of lysozyme macroions. These complexes are found to be sta- 
ble even in the presence of significant concentrations of monovalent salt, which is quantitatively 
explained from a redistribution of salt between the condensed and the aqueous phases. 

PACS numbers: 82.35.Rs, 87.16.Ka, 87.64.Bx, 87.15.Aa 
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In the presence of multivalent cations, anionic biolog- 
ical polyelectrolytes can overcome their electrostatic re- 
pulsion and exhibit a mutual attraction. These "like- 
charge attractions" result from ion correlations that can- 
not be understood within mean-field theories such as 
the commonly-employed Poisson-Boltzmann formalism 
Q, 1^ 1^ . The problem becomes more complex when the 
mediating multivalent cations are themselves macroions. 
Macroion-polyelectrolyte complexes occur in many phys- 
ical systems, such as DNA-dendrimer complexes for non- 
viral gene therapy 0] and antimicrobial binding in cys- 
tic fibrosis 0|. Various factors affect their formation: 
The presence of salt can lead to an attraction driven 
by osmotic pressure 0. Differential screening of posi- 
tive and negative charges distributed on the surface of 
a macroion may significantly modify interactions at the 
macroion-polyelectrolyte interface [7| . Entropic gain due 
to mutual neutralization and consequent counterion re- 
lease upon macroion-polyelectrolyte "adhesion" is ex- 
pected to be important, but can be potentially modu- 
lated by the steric commensurability between the charge 
pattern on the polyelectrolyte and the macroion size [Sj. 
The relative importance of all these interactions, and 
how they modify one another in their combined effect 
on the structural evolution of macroion-polyelectrolyte 
complexes, is generally unknown. 

In this Letter, we examine the role of several of 
the above-mentioned interactions in the complexation of 
actin and lysozyme, a prototypical system of oppositely- 
charged "rods" and "spheres," over a range of monova- 
lent salt concentrations. Using synchrotron small-angle 
x-ray scattering (SAXS), we show that self-assembled 
complexes are comprised of hexagonally-coordinated 
columnar arrangements of actin held together by one- 
dimensional (ID) arrays of lysozyme macroions at the 
three-fold interstitial "tunnels" of the columnar actin 
sublattice (Fig. Molecular dynamics (MD) simu- 

lations using a realistic model of the actin helix pro- 
vide a detailed confirmation of this picture, and reveal 




FIG. 1: (A) Synchrotron 2D x-ray diffraction pattern of par- 
tially aligned actin-lysozyme bundles, formed in a solution 
containing 150mM KCl. (B) ID integrated slices along the qz 
and Qr directions with arrows marking the actin-actin close- 
packed bundling peak (1), the actin helix form factor (2), and 
the lysozyme-lysozyme correlation peak (3). (C) Proposed 
structure of actin-lysozyme composite bundles (side and end 
views): Lysozyme (orange) is close-packed in three-fold sym- 
metric sites between actin filaments (blue). 



structural reconstructions and corresponding salt redis- 
tribution within an actin-lysozyme bundle as the inter- 
actin separation is varied. Both experiment and sim- 
ulation show that the lysozyme is arranged in a close- 
packed manner, incommensurate with the actin period- 
icity. Moreover, the self-assembly of columnar actin- 
lysozyme complexes is enhanced for higher concentra- 
tions of monovalent ions. We believe that these re- 
sults can be explained by significant repartitioning of salt 



2 



between the condensed and the aqueous phases, which 
strongly modifies screening effects. 

F-actin is an anionic rod-like cytoskeletal polymer (di- 
ameter 7.5nm, charge density — e/0.25nm, persistence 
length lOfim). Lysozyme is approximately an ellipsoid 
of size 2.6nm x 2.6nm x 4.5nm with a net charge of 
+9e at neutral pH. Monomeric G-actin (MW 42,000) 
was prepared from a lyophilized powder of rabbit skele- 
tal muscle. The non-polymerizing G-actin solution con- 
tained a 5mM TRIS buffer at pH 8.0, with 0.2mM CaCla, 
0.5mM ATP, and 0.2mM DTT and 0.01% NaNg. G- 
actin (2mg/ml) was polymerized into F-actin upon the 
addition of salt (lOOmM KCl). Human plasma gel- 
solin was used to control the average F-actin length 
to ~lmm. The F-actin filaments were treated with phal- 
loidin (MW 789.2) to prevent depolymerization. Hen egg 
white lysozyme (MW 14,300) was mixed with F-actin 
in 1.5mm diameter quartz capillaries to form isoelectric 
actin-lysozyme complexes. SAXS experiments were per- 
formed both at Beamline 4-2 of the Stanford Synchrotron 
Radiation Laboratory as well as at an in-house x-ray 
source. The incident synchrotron x-rays from the 8-pole 
Wiggler were monochromatized to 8.98KeV (A=1.3806A) 
using a double-bounce Si(lll) crystal, focused using a 
cylindrical mirror. The scattered radiation was collected 
using an MAR Research charge-coupled device camera 
(pixel size 79fim). For the in-house experiments, inci- 
dent CuKq radiation (A=1.54A) from a Rigaku rotating- 
anode generator was monochromatized and focused using 
Osmic confocal multilayer optics, and scattered radiation 
was collected on a Bruker 2D wire detector (pixel size 
105/im). The 2D SAXS data from both systems are mu- 
tually consistent. 

A 2D diffraction pattern for partially aligned isoelec- 
tric F-actin-lysozyme bundles and its associated ID in- 
tegrated slices along the qz and qr directions are shown 
in Figs. nj\ and ^3. Examination of the slice along 
the equatorial (qr) direction shows a correlation peak 
at q—0.07A~^ that corresponds to close-packed compos- 
ite actin-lysozyme bundles. The inter- actin spacing of 
90A is consistent with a columnar actin lattice expanded 
by lysozyme in three-fold interstitial "tunnels," aligned 
with its long axis parallel to the actin. This inter-actin 
spacing is significantly larger than the 75A spacing for 
close-packed actin condensed with multivalent ions 0. 
No other arrangement of lysozyme and actin will repro- 
duce this diffraction pattern, given their respective sizes. 
Along the meridional (qz) direction, a weak, mosaic- 
smeared actin form factor feature at O.IISA"^ is ob- 
served, as well as a new, strong correlation peak that 
differs from expected actin form factor features (layer 
lines) in position, orientation, and relative intensity. The 
appearance of this peak at qz—0-lSOK~^ corresponds 
to an inter-lysozyme distance of 48. SA, comparable to 
the length of lysozyme along its major axis, which sug- 
gests that lysozyme is close-packed along this direction 



within the bundles. Interestingly, this lysozyme periodic- 
ity is incommensurate with the projected actin periodic- 
ity (~56A), in contrast with the behavior of divalent ions 
on actin This incommensurate arrangement permits 
charge matching between actin and lysozyme within the 
bundle, and indicates the important role of entropy gain 
from counterion release in this system. Figure ^3 shows 
schematic representations of a condensed bundle. 

In order to elucidate the underlying mechanism for 
bundle formation and the structure of the resulting com- 
plex, we have performed MD simulations using a mod- 
ified version of Moldy 0|. In these simulations, G- 
actin is modeled using the four-sphere model 0|, which 
is based upon crystallographic measurements and pro- 
vides a relatively accurate coarse-grained representation 
of the monomer charge distribution. F-actin is comprised 
of a sequence of these monomers, in which successive 
units have a separation of 27. 5A and a relative rota- 
tion of 166.7° around the filament axis. This leads to 
a helical structure with a repeat unit of 13 monomers. 
The filaments are assembled into a parallel hexagonally- 
coordinated bundle. An elementary simulation cell con- 
sists of a bundle fragment containing 2x2 filaments with 
a length of 6 repeat units (78 monomers) each. This cell 
is periodically replicated in all directions. Following the 
experiments, we set the lysozyme concentration in the 
bundle to neutralizing conditions, which corresponds to 
352 lysozyme units per simulation cell. Each lysozyme 
is modeled as a rigid dumbbell structure of two spheres 
with diameter 25A and charge 4.5e, at a center-to-center 
distance of 20A, thus approximating the aforementioned 
ellipsoidal dimensions. Additional salt is modeled as 
monovalent spherical particles with a hydrated radius 
of 3.3A. Coulomb interactions are treated by means of 
Ewald summation, and excluded-volume interactions are 
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FIG. 2: Osmotic pressure of a hexagonally-coordinated bun- 
dle of actin filaments in excess solution, as determined from 
MD simulations. Inset: lysozyme pair correlation function 
along the filament axis. For discussion see the text. 



represented by pair potentials of the form k^T{a /r)^"^ , 
where a is the sum of the effective radii of two interact- 
ing particles (ions, G-actin subunits and lysozyme sub- 
units). During each simulation, the actin separation is 
fixed, whereas lysozyme and all ions move freely. Thus, 
the simulations probe the stability of a swelling bundle 
while maintaining the filaments in a parallel configura- 
tion, ignoring their rotational degrees of freedom. This 
is justified by the observation that the calculations are 
confined to actin separations below 25nm, i.e., less than 
0.25% of the persistence length. Mutual sliding and axial 
rotation of the filaments are not taken into account. 

This model indeed predicts electrostatically driven 
complex formation. Since the water is modeled as a di- 
electric continuum, the osmotic pressure 11 can be ob- 
tained directly from the virial involving all interparticle 
forces Bundle formation takes place in excess solu- 
tion, and hence the bundle stability follows from a com- 
parison of the osmotic pressure to the osmotic pressure 
of the salt Ilsait- A negative osmotic pressure difference 
An = n — Ilsait implies bundle contraction and the free- 
energy minimum [All = and d{/^I\)/dV < 0] yields 
the stable actin separation. A coinparable approach has 
been employed before (see Refs. 0, and references 
therein) to study the condensation of rod-like polyelec- 
trolytes by counterions. FigureElshows that, under salt- 
free conditions, an inter-actin spacing of '^lOOA is pre- 
dicted, in quite close agreement with the experimental 
observations. 

Having established that our model captures essential 
aspects of actin-lysozyme complexation, we exploit it 
to elucidate the structural properties and evolution of 
the resulting complexes. Figure |21 shows contour plots 
of the lysozyme center of mass, projected on a plane 
perpendicular to the bundle axis. In order to mini- 
mize artifacts resulting from the finite actin length, the 
calculations employed filaments consisting of 12 repeat 
units (156 monomers). In the equilibrium configuration 
(Fig. the maximum lysozyme concentration occurs 
in the three-fold interstitial regions between the actin fila- 
ments, supporting the interpretation of the SAXS diffrac- 
tion data in Fig. Furthermore, the lysozyme pair cor- 
relation function (inset of Fig. along the actin axis 
shows a clear peak at z — 50A, corroborating the ex- 
perimentally measured close-packed value of 48. 3A [l5| . 
Minor enhancements of the lysozyme concentration can 
also be seen in the bridging regions between pairs of 
neighboring filaments. The lysozyme distribution is gov- 
erned both by entropic effects and by a competition be- 
tween electrostatic protein repulsions and actin-lysozyme 
attractions. Also, the helical actin structure imposes 
an excluded-volume repulsion with a rather pronounced 
shape. In order to distinguish any structure resulting 
from this repulsion. Fig. |33 shows a contour plot ob- 
tained in a simulation without electrostatic effects. The 
homogeneous lysozyme distribution confirms that the 
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FIG. 3: Contour plots showing the lysozyme distribution in 
actin-lysozyme complexes without added salt. Darker shad- 
ing corresponds to higher concentrations. In the equilibrium 
configuration (A), lysozyme is predominantly located in the 
three-fold interstitial regions. The free-energy inflection point 
occurs at a slightly expanded lattice (C), in which lysozyme is 
depleted from these regions and occupies the bridging regions 
instead. Panels (B) and (D) are the counterparts of panels (A) 
and (C), respectively, in the absence of electrostatic forces. 



structure in Fig. I^JV is dominated by electrostatic inter- 
actions. Interestingly, the actin excluded volume (white 
circular regions) also has a more pronounced hexagonal 
structure in panel A, which is caused by the precessing 
highly-charged regions on the monomers. Consideration 
of larger lattice spacings provides information on inter- 
mediate states that may arise during the complexation 
process. As illustrated in Fig. significant rearrange- 
ments occur in the final stages of the bundle formation: 
At the osmotic-pressure minimum (free-energy inflection 
point), i.e., at an actin separation that is increased by 
merely 16A, lysozyme is depleted from the interstitial re- 
gions and instead predominantly occupies the bridging 
sites. Figure|3jl) confirms that this lysozyme distribution 
is again dominated by electrostatic effects. Thus, the 
complex evolves from lysozyme at bridging positions be- 
tween pairs of actin rods to the final three-fold positions 
observed by experiment. Lysozyme and actin "contact" 
interactions are maintained during this structural relax- 
ation, indicating that their mutual electrostatic attrac- 
tion plays an important role. 

It is interesting to consider how this self-assembly is 
affected by the addition of monovalent salt. The counter- 
ion release mechanism implicit in actin-lysozyme bind- 
ing will strongly modify qualitative arguments based on 
screening. The repartitioning of salt between the con- 
densed and aqueous phases implies a different degree of 
screening inside and outside the complex. In addition, it 
can lead to a stabilizing external osmotic pressure. Such 
a redistribution of ions is often not taken into account 
theoretically. To study these effects, a series of SAXS 
measurements were performed on actin-lysozyme com- 
plexes at different NaCl concentrations. As [NaCl] is 
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FIG. 4: Series of difFraction data showing the evolution of 
bundle structure as a function of (A) NaCl and (B) KCl con- 
centration with maximum bundling occurring around 120mM. 
(C) Simulated (o) bundle vs. bulk ion concentrations for a sta- 
bilized actin-lysozyme complex. Salt repartitions into differ- 
ent concentrations inside and outside the bundle. The dashed 
line is a guide to the eye. 

increased to ISOmM, the turbidity increases and the in- 
tensity of the lysozyme-actin difFraction peak increases 
without significant changes in its peak width (Fig. 
indicating the formation of more bundles (rather than 
bundles that are more ordered and have larger coherent 
domains). At higher salt concentrations, the trend re- 
verses and a weakening of the bundling peak is observed. 
The same results are found using KCl (Fig. 0)3), showing 
that this is not a cation-specific binding effect. 

Since regular screening is likely to play a role in the 
ultimate disappearance of the bundle, we concentrate on 
the stability at low and intermediate salt levels. Simu- 
lations of bundles with additional salt show that the os- 
motic pressure within the complex rises more rapidly as 
a function of salt concentration than the bulk pressure of 
salt at the same concentration, leading to destabilization 
once the salt concentration exceeds ~10mM. This appar- 
ent discrepancy with the experimental findings already 
suggests a redistribution of salt ions. To quantify this 
further, grand-canonical simulations of the bundle can be 
employed Here, we have chosen an alternative strat- 
egy. Using the Widom particle-insertion technique (l6l |. 
we determine the chemical potential of salt within the 
complex, as a function of concentration and actin sep- 
aration. Exploiting the coexistence condition, we sub- 
sequently determine the corresponding bulk salt concen- 
trations and osmotic pressures via independent grand- 
canonical Monte Carlo simulations. This yields several 
important results. Firstly, the salt concentration inside 
the actin-lysozyme complex is approximately twice lower 
than the bulk concentration, which leads to a difference 
in osmotic pressure that is sufficient to maintain bun- 
dle stability up to much higher concentrations than oth- 
erwise would have been possible (Fig. ^p). Secondly, 
the depression of the ion concentration within the bun- 



dle may explain why maximal actin-lysozyme bundling 
is observed at a global salt concentration around 120mM, 
rather than at concentrations in the range 30-60mM as 
predicted from differential screening arguments "Tl . 

In summary, by studying self-assembled actin- 
lysozyme complexes via a combination of small-angle x- 
ray scattering and molecular dynamics simulations, we 
have shown that salt repartitioning impinges strongly on 
the structure and stability of the complex, and quali- 
fies commonly-invoked mechanisms such as counterion 
release and differential screening. 
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